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Introduction: Before a hydrocode is used to 
investigate a question of scientific interest, it should 
be tested against analogous laboratory experiments 
and problems with analytical solutions. The 
Radiation Adaptive Grid Eulerian (RAGE) 
hydrocode[l], developed by Los Alamos National 
Laboratory (LANL) and Science Applications 
International Corporation (SAIC)[2,3] has been 
subjected to many tests during its development. [4,5] 
We extend and review this work, emphasizing tests 
relevant to impact cratering into volatile-rich targets. 

RAGE: is an Eulerian radiation-hydrodynamics 
code that runs in a variety of geometries in up to 
three dimensions, with a variety of equations of state. 
It was developed for general application, so it does 
not include ad hoc tuning of algorithms or 
parameters, but relies solely on physical first- 
principles. RAGE uses a higher order piecewise 
linear Godunov numerical method to solve the 
hydrodynamics equations. [6,7] An exceptional 
advantage to this code is its use of a time- and space- 
continuous adaptive mesh refinement (AMR), by 
which it is able to follow shocks and other 
discontinuities at high resolution while treating 
smooth regions coarsely, increasing computational 
efficiency. The radiation component of the code is 
an optional grey diffusion model with nonequilibrium 
radiation and material temperatures. [8] Opacities are 
generated as part of the SESAME equation of state 
tables, or optionally from an analytic model. 

Verification. This process ensures that the model 
is coded and solved correctly. RAGE is subject to an 
extensive set of tests in which output of each version 
of the code, on every computer system on which it is 
run, is compared to analytical solutions and results of 
previous versions of the code in order to demonstrate 
invariance under different running and boundary 
conditions. Test problems include the Sod shock 
tube[9] demonstrated in this study, the Noh problem, 
which checks for errors from shock 'smearing' by 
finite resolution] 10], the Sedov blast wave, which 
scales self-similarly in timefll], Marshak waves, 
which tests radiation diffusion [12], and many other 
problems [4], 

Validation. This process ensures that the model is 
appropriate to the problem at hand, and of sufficient 
accuracy. To validate the RAGE code, many 
simulations are used to replicate physical 


experiments. These include examinations of fluid 
instability in a shock-accelerated thin gas layer[13], 
Richtmyer-Meshkov instability growth[14], 
supersonic fluid flow and shock-induced jetting[ 15], 
and shock transmission through boundaries [16]. We 
extend this effort to shocks in basalt and ice. 

SESAME: is a temperature-based tabular 
equation of state maintained by the Mechanics of 
Materials and Equations of State group LANL. The 
table for each material has a unique and 
thermodynamically consistent fit of semi-empirical 
theoretical models appropriate to different 
temperature or pressure regions to experimental 
datafl 7, 1 8]. 

Method: We conduct four simulations of 
particular relevance to impact modeling. The first is a 
reproduction of the one-dimensional shock tube 
verification problem published by Sod[9], Initially 
the 1-D planar tube is divided into two sections. In 
the left section, a gas is in equilibrium at a higher 
density and pressure (pi=1.0, pi=1.0, Ui=0.0). In the 
right section, the gas is in equilibrium at a lower 
density and pressure (p:=0.125, p:=0.1, u:=0.0). At 
t=0 the high-pressure/density gas expands into the 
low-pressure/density region, generating a shock 
wave. The hydrocode results are compared to the 
analytical solution to demonstrate the accuracy of the 
finite differencing scheme. The second and third tests 
reproduce shocks in basalt[19] and water ice[20]. 

Results: The Sod shock tube performed well. 
Further tests are ongoing. Detailed results will be 
presented at the workshop. 

Sod Ideal Gas Shock Tube. This model 
corresponds to the shock tube described above[9], 
and compared to the analytical solution. Results were 
graphed for density, velocity, pressure, and specific 
internal energy at t=0.12 seconds (fig. 1-4), just 
before the shock wave hit the end of the tube. 

Ice Shock Tube. In this experiment[20], aluminum 
and polycarbonate projectiles were fired into water 
ice slabs of thicknesses ranging from 3-60 mm. They 
measured pressure, attenuation, and propagation 
velocity. In RAGE, we will use a one-dimensional 
column where an aluminum[21] or polycarbonate[22] 
projectile strikes an ice (SAIC proprietary water 
equation of state) target at t=0, at velocities ranging 
from 300-600 m/s, as in Kato[20], Pressures and 
velocities will be recorded for comparison with the 
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original data. 

Basalt Shock Tube. In this experiment 19], a 
copper projectile was fired at a column of basalt 
plates interleaved with pressure sensors, generating 
peak shock pressures from 7-9 GPa. They measured 
the pressure over time through the column. In RAGE, 
we will use a one-dimensional column where a 
copper[23] projectile will strike a basalt[24] target at 
t=0, at velocities between .68-2.7 km/s. Shock 
pressure and propagation will be recorded for 
comparison with the original data. 

Future Work: When we are confident that 
verification and validation criteria have been met, we 
will use the RAGE hydrocode to examine in two 
dimensions the effects of sub-surface volatiles, such 
as water ice, on the propagation of impact-generated 
shock waves through the Martian surface. 

References: [1] Gittings, M. L., et al., (2004) 
computer code RAGE ver. 20041126.002 (LANL, 
SAIC). [2] Clover, M. and Gittings, M. L. (2002) 
APS DCOMP 2002, G2001. [3] Clover, M. (2002), 
LA-UR-02-3801. [4] Kamm, J. R. and Rider, W. J. 
(2000), LA-UR-99-5234. [5] Clover, M. (1998), LA- 
UR-98-5570. [6] Woodward, P. and Colella, P. A. 
(1984) J. Comp. Phys, 54, 115-173. [7] Miller, G. H. 
and Colella, P. A. (2001), J. Comp. Phys, 167, 131- 
176. [8] Saatdjian, E. (2000), Transport Phenomena: 
Equations and Numerical Solutions , pp. 181-186, 
(Wiley, New York). [9] Sod, G. (1978), J. Comp. 
Phys, 27, 1-31. [10] Noh, W. F. (1978), J. Comp. 
Phys, 72, 78-120. [11] Castor, J.I. (2004), Radiation 
Hydrodynamic s, pp. 35-37, (Cambridge). [12] 
Zel'dovich, Ya. B. and Raizer, Yu. P. (1967) Physics 
of Shock Waves and High-Temperature 
Hvdrodvnamic Phenomena , v. 2, pp. 672-676, (AP, 
New York). [13] Baltrusaitis, R. M. et al. (1996), 
Phys. Fluids, 8(9), pp. 2471-2483. [14] Holmes, R. L. 
et al. (1999), J. Fluid Mech., 389, pp. 55-79. [15] 
Foster, J. M. et al., (2002), Phys. Plasmas, 9(2), pp. 
2251-2263. [16] Goldman, S. R. et al. (1999), Phys. 
Plasmas, 6(8), pp. 3327-3336. [17] Johnson, J. D. 
(1994) 1 2th Symp. Thermophysical Properties. [18] 
Lyons, S. P., and Johnson, J. D. (1992), LA-UR-92- 
3407. [19] Nakazawa, S. et al. (2002), Icarus, 156, 
pp. 539-550. [20] Kato, M. et al. (2001), JGR, 
106(E8), pp. 17,567-17,578. [21] Chisholm, E. D. et 
al. (2002), LA-UR-6691. [22] Boettger, J. C. (1991), 
LA-12 120-MS. [23] Hayes, D. B. et al. (1999), LA- 
UR-3276. [24] Bames, J. F., and Lyon, S. P. (1988), 
LA-1 1253-MS. [25] 

Acknowledgments: UCSC/LANL Institute of 
Geophysics and Planetary' Physics, NA.SA PG&G 
Small Bodies and Planetary Collisions. 



Figure 1. Sod ideal gas shock tube density vs. 
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Figure 4. Specific internal energy vs. distance, 
t=0.15s. 






